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In an attempt to derive thermodynamics from classical mechanics, an approximate expression for 
the equilibrium temperature of a finite system has been derived [M. Bianucci, R. Mannella, B. J. 
West, and P. Grigolini, Phys. Rev. E 51, 3002 (1995)] which differs from the one that follows from 
the Boltzmann principle S = k\nQ(E) via the thermodynamic relation 1/T = dS/dE by additional 
terms of "dynamical" character, which are argued to correct and generalize the Boltzmann principle 
for small systems (here Q(E) is the area of the constant-energy surface). In the present work, the 
underlying definition of temperature in the Fokker-Planck formalism of Bianucci et al. is investigated 
and shown to coincide with an approximate form of the equipartition temperature. Its exact form, 
however, is strictly related to the "volume" entropy S = k In <&(E) via the thermodynamic relation 
above for systems of any number of degrees of freedom (<fr(E) is the phase space volume enclosed 
by the constant-energy surface). This observation explains and clarifies the numerical results of 
Bianucci et al. and shows that a dynamical correction for either the temperature or the entropy is 
unnecessary, at least within the class of systems considered by those authors. Explicit analytical 
and numerical results for a particle coupled to a small chain (TV ~ 10) of quartic oscillators are also 
provided to further illustrate these facts. 
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I. INTRODUCTION 



In the notation of Planck, the connection between the thermodynamic entropy S of an isolated system and its 
number of accessible states W is often introduced as what Einstein called the Boltzmann principle pj 

S = k\nW. (1) 

For classical systems, where the number of states is replaced by continuum integrals over the phase space, two common 
definitions of entropy for isolated systems are 1 (see e.g. @, or Q and references therein) 

S$ = k In $(#), (2) 
Sn = klnQ{E), (3) 



where 



is the volume under the constant-energy surface, and 



$(£)= \ dz = [ dz6[E - H(z)] 



(4) 



Sl(E)= [ dz= [ dzS[E- H(z)} (5) 

is the "area" of that surface. Here H(z) is the Hamiltonian, z denotes a point in the phase space, 0{x) is the 
Heaviside step function and 8(x) is the Dirac delta function. These two expressions for S will be referred to as the 
volume and area entropies, respectively. It is generally acknowledged that a specific choice of S above is immaterial 
in the thermodynamic limit 0j 13 > for therein these expressions differ by terms increasingly small, but in the opposite 
limit of only a few degrees of freedom they lead to appreciably different observables, most notably the equilibrium 
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1 The prefactors 1/7V! and l/h' iN are omitted since they are irrelevant in the present context. 
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temperature (given by the thermodynamic relation 1/T = dS/dE) Q. The temperatures that follow from Eqs. (j2J 
and via this thermodynamic relation will be similarly referred to as the volume (T$) and area (Tq) temperatures. 

Almost a decade ago, Bianucci et al. (BMWG henceforth) set out to provide a "derivation," and generalization, 
of the Boltzmann principle 2 [Eq. @] with particular emphasis on systems far away from the thermodynamic limit 
Such generalization, which will be concisely reviewed in the next section, is based upon the derivation of an 
approximate equation of motion for the reduced probability distribution of the degrees of freedom of interest, which 
is then compared with a postulated form of the Fokker-Planck equation [Eq. i|12fl below]. This comparison allows 
BMWG to put forward a microscopic expression for the temperature [Eq. I|14|) in the present work], and additional 
assumptions allow them to relate their temperature expression Tb with the phase space area Q [Eq. I|18p]. Since the 
resulting relationship between Tg and f2 differs from the one that follows from Eq. @ via the thermodynamic relation 
1/T = dS/dE by additional terms that only vanish in the thermodynamic limit, Bianucci et al. are led to conclude 
that they have found a generalization of the Boltzmann principle for finite systems. This fact is confirmed by BMWG 
through numerical simulations of a particle coupled to a chain of oscillators, in which the authors measure the average 
of twice the kinetic energy of the particle and compare it with the predicted value from their proposed temperature 
expression. 

It is clear from the above discussion that an implicit assumption in the work of BMWG is the coincidence of their 
analytic temperature definition Tb [Eq. l|18f) ] with (v 2 ') (twice the average kinetic energy, where we temporarily set 
the mass and the Boltzmann constant to unity), i.e. it is assumed that their temperature definition coincides with, 
or at least approximates the familiar equipartition one. Since it was actually the latter quantity that was used to 
assess the quality of the former and not the other way around, one can say that the endeavor of BMWG was to find 
a dynamical method for evaluating (v 2 ) in terms of the phase space structure of the system, and that the failure of 
Eq. J2J in predicting the equilibrium temperature of their numerical experiments is related to the breakdown of the 
equation (u 2 ) = (dSn/dE)" 1 for finite systems. As shown in the present investigation, the BMWG temperature Tb is 
in fact an approximation to the equipartition one, but if one assumes mixing as BMWG did in the final stages of their 
investigation, the asymptotic value of (u 2 ) becomes exactly related to S$ [Eq. J21J] for any system size, establishing 
incidentally the equivalence between the equipartition and the volume temperatures. This shows not only that the 
Boltzmann principle in Eq. @ is indeed inappropriate to predict the equipartition temperature of small systems Q , 
but also that the dynamical corrections to So suggested in the work of Ref. [j| are unnecessary, since an alternate and 
equally legitimate form of the Boltzmann principle [Eq. @] is sufficient for that purpose. This observation translates 
immediately into a negative answer to the question posed in the title of the present work, so long as one answers it 
in the context of the work of BMWG. 

The remaining of this paper is organized as follows. After reviewing the formalism of BMWG (Sec.[n]), it will be 
shown that the temperature expression found in Ref. p| is an approximation to the usual one defined through the 
equipartition theorem, and hence that the expression derived therein as a generalization of the Boltzmann principle 
is but a particular approximation to the volume entropy fSec. These findings will be illustrated by means of a 
Hamiltonian model similar to that of BMWG (Sec. HV() . where in addition the distribution of momenta is computed 
and seen to depart from the canonical one. This provides, in turn, an alternate means of explaining the inability of the 
area temperature to predict the average of twice the kinetic energy of a particle, as seen in the numerical experiments 
of Bianucci et al. 



II. THE DYNAMICAL APPROACH OF BMWG 



In order to understand the origin of the dynamical generalization of the Boltzmann principle mentioned above, it 
is worth reviewing the scope and results of BMWG. The authors of Ref. [f| were concerned with problems whose 
degrees of freedom can be separated into "system of interest" and "irrelevant," the latter playing the role of what 
is commonly known as heat bath or thermostat in the literature. However, in order to avoid confusion with the 
usual Fokker-Planck equation approach to their problem, in which the irrelevant degrees of freedom are assumed to 
be in a specific macroscopic state (usually a canonical one, see e.g. 0,0), the word booster was adopted instead. 
Therefore, the general goal of BMWG was to derive the thermostatistical properties of a system coupled to a (finite 
or infinite) bath whose thermal properties arise naturally from dynamics, bridging the gap between mechanics and 
thermodynamics much like the ideas of Boltzmann. 

Notwithstanding the avoidance of macroscopic assumptions, two fundamental properties concerning the dynamics 
of the booster had to be assumed in order to make the problem amenable to their analytical approach, namely (i) it 
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is chaotic enough so that its correlation functions decay in a finite time (finite correlation time) and (ii) the average 
response of the system to an external perturbation is linear (linear response). In addition, the system of interest was 
taken to be a single particle with coordinates x, v linearly coupled to the booster via a "doorway" variable £, so that 
A • £ can be interpreted as an external force driving the particle (or, vice versa, A • x for the booster) whose coupling 
strength A has to be sufficiently small to justify their perturbative approach (weak coupling condition). Condition (i) 
implies in particular that the unperturbed autocorrelation function of the doorway variable (with (-) indicating an 
equilibrium average in the absence of coupling), viz. 

**) = -g^, (6) 

vanishes for finite t (with (£} = 0), which in turn implies that the correlation time 

poo 

dt ip(t) (7) 



f 

Jo 



is finite. Condition (ii) states that the average behavior of the doorway variable after the introduction of an external 
time-dependent perturbation K(t) to its equation of motion can be written as the following convolution (recall the 
assumption (£) = above) 



t 

2 



<£(*)>*= / du S(u)K(t-u) + 0(K z ), (8) 



where S(u) is a function that dictates how the average of £ responds to the small perturbation K(t) (response 
function). Together with (i), condition (ii) furnishes two other finite quantities, the asymptotic susceptibility 



t 

du S(u), (9) 



X = lim x(t) = lim / 

t — >oo t^-oo Jq 

and the response time 

nOC 

§= I dtc(t), (10) 



where the convenient variable c(t) = 1 — was introduced so that c(0) = 1 and c(oo) = 0. 

Relying on the above assumptions (and also on a few others concerning the separation of time-scales) and equipped 
with a specific projection technique, Bianucci and his collaborators embarked on the rather remarkable endeavor of 
"integrating out" the irrelevant degrees of freedom of the booster without the precise knowledge of its Hamiltonian 
structure, ultimately in the search of a differential equation governing the time evolution of the reduced distribution 
function of the particle, which can be formally written as |5( 

cr(x,v;t) = J d^d-K p(x, v, n, £; i), (11) 

where p(x, v, 7r, £; t) is the distribution function of all the dynamical variables and tt denotes the remaining degrees 
of freedom of the booster. The resulting differential equation, Eq. (47) in Ref. |5(, is then identified with a general 
2D Fokker-Planck equation of the type 

d -a(x,v;t) = [cf+^-A(x,v)(- d 



dt v ' ' ' \ a dv y ' ' \dv kTj 
d J d U'(t 



+ ±B(x,v)(^-+ l ^)}a(x,v;t), (12) 



where 



£f=^|--4 (13) 
m av ox 

is the effective time evolution operator (Liouvillian) of the unperturbed particle. From a comparison between the 
coefficients of their equation and the structure of the above equation, they suggest (despite the implicit dependence 
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on the particle coordinates x, v) that the following expression should be identified with the temperature of the booster 
(Eq. (52) in Ref. @): 

k T B (x,v) - ^J U ^ {u) ^ Xait ~ u) ] , (14) 

X f™duc(u) [£a: (t-u)] ' 

where x a (t — u) is the unperturbed evolution of the position of the particle "backwards in time" (being essentially 
the responsible for the implicit dependence of Tb on x,v). The subscript in Tb stands for BMWG. To render the 
above equation independent of x, v, BMWG consider the following three cases: (a) the particle is a harmonic oscillator 
with V(x) = mui 2 x 2 /2, (b) (p(u) = c(u), and (c) a "natural time scale" (NTS) condition exists separating the typical 
timescale of the system of interest from the relaxation TS of the booster (the former being much greater than the 
latter). Under condition (a), the temperature reduces to 

_ <a«a, (15) 

with hats indicating a Fourier transform, whereas in the case of NTS (c) one obtains 

kT B = ^4-' (16) 

which coincides with the one obtained within a Langevin approach (cf. Eq. (22) in Ref. @), as expected from the 
dichotomy of the time scales. Condition (b), however, is peculiar in the sense that it recovers the usual canonical 
results, which (with the exception of making the connection with Kubo's linear response theory clear) was not 
particularly relevant for the remaining part of their work. 

To draw a parallel with Boltzmann's principle, BMWG invoked mixing for the booster, in which case the time- 
dependent susceptibility can be found via "geometrical" arguments to assume the form H 

^ = ^m^^^-m}, d7) 

where Slb(E) is the phase space area of the unperturbed booster with energy E. With this result in hand, one can 
write, for the case (c) above for example, 

^ = ^ 1 " n >< £ > + l? 1 "««V>. us) 

whereas for case (a) one simply replaces r with Re [^(w)] in this expression. This result should be contrasted with 
the one obtained from what BMWG calls the Boltzmann principle [i.e. Eq. J2J], namely 

lk = JE^ E) - (19) 

The additional term in Eq. (|18(l . called by BMWG a "dynamical correction" to Eq. (|19|l . is then argued to become 
negligible for sufficiently large systems, thus recovering the result above and causing no conflict with the standard 
macroscopic thermostatistics of Boltzmann, i.e. the Boltzmann principle has been "generalized" to cover systems with 
only a few degrees of freedom. In the following section it is shown that this generalization is actually an approximate 
form of the volume temperature T$ . 



III. CONNECTION WITH THE VOLUME TEMPERATURE AND ENTROPY 

Besides heuristic arguments of units, it is not clear how one can identify the "constant" Tb [Eq. I)14[l] obtained by 
BMWG through a comparison with Eq. I|12(l as an actual temperature. In the usual derivation of the Fokker-Planck 
equation (or, equivalently, of the Langevin equation) for a Brownian particle 0, ITo| . the temperature comes into 
play by invoking the equipartition theorem, viz. kT = m(v 2 ) (this is indeed the case for the Langevin treatment of 
BMWG, cf. Eq. (22) in Ref. @), but in the case of the Fokker-Planck equation derived by BMWG or the generic 
one in Eq. (|12|l . both with x, ^-dependent coefficients, it is difficult to see whether their temperature definition is 
compatible with equipartition. Nevertheless, after the simplifications introduced by any of the cases (a), (b), or (c) 
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of the previous section (the only ones considered by BMWG in their generalization of the Boltzmann principle) , one 
can show that the definition of kTg coincides with an approximate form ofm(v 2 ). Indeed, focusing on condition (c) 
for notational simplicity (the remaining cases (a) and (b) follow by a trivial modification of the procedure below), 
the time-independent differential equation satisfied by a(x,v) obtained via the projection method reduces to (cf. Eq. 
(47) of Ref. H) 

da V'(x)da A 2 x X da A 2 (£ 2 ) r «9 2 cr 
= — v 1 — — 1 — — 

dx m dv m dv m 2 dv 2 

A 2 (£ 2 ) 7? 2 d 2 ~a AV d 



m 2 dvdx m dv 
A 2 x /3 2 da 
m dv ' 



(va) 

(20) 



where the constants rj 2 and 1 are defined in Eqs. (6) and (13) of Ref. p|- The tilde indicates that the distribution here 
is not identical to Eq. due to the approximations carried out in the actual process of integration, or "projection" 
(see below). Multiplying the above equation through by v 2 and integrating over x, v, most terms disappear under the 
assumption that a(x,v) — » sufficiently fast for x,v — > oo and (x)& — {v)& — (for simplicity I assume the average 
position of the particle is centered at the origin, though it is not difficult to adapt this derivation for the general case), 
so one is left with 

= ^ ^(v 2 ) a , (21) 

where (■)& = dx dv ■ a(x, v). Solving this equation for m (v 2 ) ~ we obtain the right hand side of Eq. (TSJ, 
i.e. m(v 2 ^^ — IzTb- This is the desired connection between Tg and the equipartition temperature, q.e.d. 

It is important to emphasize, however, that the reduced density a (and hence the aforementioned connection) is 
only approximate. This can be seen in the discussion of paragraph 1, page 3005 of Ref. |{| - apart from conditions (i) 
and (ii) of the previous section, it is always assumed that there is a large separation of relaxation and intrinsic time 
scales of the booster. This condition, which amounts to an assumption of weak coupling, was used explicitly between 
Eqs. (39) and (40) of Ref. [f| (note also that Eq. (39) of Ref. is itself a perturbative result). We now note that the 
exact average m(v 2 ) a , obtained with a [Eq. (|ll|l ] before the approximations carried out in the projection method of 
BMWG which lead to a, furnishes the following expression 

2 X _ L j 2 6[E-H(x,v,£,Try 



m{v 2 ) a = m / dxdvd^dir v 2 -± )' n , (22) 

J 11(E) 

where H(x,v,^,tv) is the Hamiltonian of the system + booster, and 

5[E - H(x,v,£,tt)] 
P^ v ^ = U(E) (23) 

is the asymptotic solution of the Liouville equation for mixingsystems (microcanonical distribution, cf. Q), in 
accordance with the assumption of mixing in Sec. IV-D of Ref. |5j. By means of simple integral manipulations used 
in the generalized equipartition theorem (see e.g. 0j0|), we can rewrite Eq. (|22|l in terms of the volume entropy 
5$ = fcln$(£0 as 



^(E) 



k(^±) =A7V 124) 



Observe that, contrary to the approximate conclusion made after Eq. 1)21(1 . the above equation establishes an exact 
equivalence between the equipartition and the volume temperatures, as already anticipated in Sec^ showing that the 
knowledge of the phase space volume alone (or S$) is sufficient to predict the equilibrium value of m (f 2 ) cr - Though 
not new 0, Q , this result renders the dynamical corrections introduced in Ref. unnecessary. This fact is illustrated 
in the next section by means of explicit analytical and numerical calculations. 



IV. QUARTIC OSCILLATORS: AN ANALYTICALLY TRACTABLE MODEL 



In this last section a specific Hamiltonian model will be analyzed both analytically and numerically. Since BMWG 
relied on a similar model to illustrate their findings, the present analysis serves the similar purpose of illustrating on 
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an equal footing the results established in the last section. It also shows that, for such small systems where Eqs. 
and (j3J) yield appreciably different results, a noncanonical distribution of momenta is a necessary and fundamental 
property, in contrast to what is argued in Ref. 0- 

The model adopted herein is described by the Hamiltonian structure 



H{z) = H x {x,p) +# 2 (0,7r) +gH lt2 (x, ^ ), 



(25) 



where x,p and <p = (0o, ■ ■ ■ , 4>n-i), t — (^o, ■ ■ ■ , T/v-i) are the particle and booster canonically conjugate coordinates, 
respectively, g is a controllable coupling constant and 



P 2 1 2 2 
Hi{x,p) = — + - LV X , 

N-l 2 N ~ 1 r^4 

1=0 i=l 



\41 



(26) 

(27) 
(28) 



are the particle, booster and coupling Hamiltonians, respectively. The intrinsic time scale of the particle is set by 
uj^ 1 , whereas its coupling with the booster is set linearly through the "doorway" variable <pQ. It should be noted that 
this model can be rendered unphysical when to — > (i.e. in the free particle limit) since the energy would no longer be 
bounded from below (due to the remaining linear coupling term Hi >2 ) and thence no stable thermodynamical state 
would exist (that shall not concern us here since ui will always be non-zero). 

The essential difference between the above model and the one adopted by BMWG is the absence of harmonic terms 
in the booster Hamiltonian, which permits a simple analytic treatment based on scaling arguments, as we shall see 
in the next sections. While not particularly relevant for the present analysis, it will also be shown that this model 
retains most of the quantitative features of the model of Ref. |5j] (see, in particular, q and £ below). More important 
than these quantitative similarities is its mixing behavior which, just as in the model of BMWG, is suggested by the 
numerical experiments through the observed finite relaxation time towards a unique equilibrium state. 

For the numerical simulations to be reported soon, the equations of motion were evolved with an optimized fourth- 
order symplectic integrator [l2| . resulting in numerical energy fluctuations AE/E smaller than one part in 10 6 . Initial 
conditions were picked randomly from the density p(x,p, <p, tv; t = 0) = 8(x)8(jp)8{H 2 — E) p|, and no dependence 
of the observables in equilibrium on these configurations was observed (being already a good indication of ergodicity 
by Birkhoff's theorem jij). The typical integration time used was t = 10 8 , after allowing a suitable relaxation time 
of order t re i = 10 7 . Since we are not interested in nonequilibrium features, no ensemble average was done (the 
coincidence of time and phase space averages relying on the Birkhoff-Khinchin theorem for ergodic systems Q). For 
simplicity, the harmonic frequency of the particle, ui, was set to unity (in contrast to the approach of BMWG, the 
specific value of u> is immaterial for the validity of the method adopted here) . 



A. Momentum distribution 



As a first probe of ergodicity, one can test whether the observed histogram of the momentum of the particle f(p) 
approximates the distribution law obtained from the structure function Q p (E) of the remaining degrees of freedom x, 
4> and 7r (see below), i.e. for large enough integration times one expects 

» p (g-p 2 /2) 

m -* m ■ (29) 

where fl(E) is the structure function of the whole system. Most quantities of interest here, including fl p , can be 
obtained from the (2A^ + l)-dimensional phase space volume 

%(E) = Jdxd(j)dTv 6(E- [x 2 /2 + H 2 (4>,tt)}) , (30) 

where 9 is the unit step function and the coupling Hamiltonian H\ 2 was neglected (weak coupling approximation). 
In fact, assuming the above volume to be bounded and sufficiently smooth for all E > (such assumption will be 
adopted throughout this section), the phase space area £l p can be obtained from <E> p via the general relation 01 



(31) 
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FIG. 1: Logarithmic plot of the histograms f(p) obtained from the numerical simulations for N = 4, 8 and g — 0.1. The slope 
of each curve is indicated and should be compared with the values predicted by Eq. 1341 . viz. 2.5 (N = 4) and 5.5 (N — 8). 
The finite-time sampling manifests here through the departure from the predicted power-law at the tail of the distributions, 
which nevertheless tends to smooth out as t — > oo. 



where i labels the desired (sub)volume. Though evaluating Eq. (|30[l can be a hard task for general Hamiltonians, for 
the choice presented in Eqs. (|26[) - l|28|l a simple scaling argument can be used to show that (cf. Ref. 01 ) : 

%(E) = ~aE^ N+2 V\ (32) 

where henceforth the symbols a, b, etc. will refer to any function that does not depend on E. Equation (|31[) then 
gives 

n p (E) = &£( 3Ar - 2 )/ 4 , (33) 
and therefore, by inserting this relation in Eq. (|29[). one gets the desired distribution law: 

g2X (3JV-2)/4 

Z y E2 

where Z is a normalization constant. This expression is valid for any N > 0, provided the dynamics of the system 
is sufficiently ergodic (recall that f(j>) was defined as a histogram). Note also that in the limit N — > oo with E/N 
constant one recovers the usual Boltzmann factor exp(— f3p 2 /2), with (3 = 2>N /AE. 

Figure ^ shows the momentum distribution function of the harmonic oscillator obtained from the numerical sim- 
ulations for the booster model in Eq. (|27|l . from which one obtains good evidence for ergodicity given its clear 
agreement with Eq. J23J. The "departure from Gaussianity" introduced by BMWG was also computed, namely 
q = 1 — (p 4 )/[3(p 2 ) 2 ], to show not only that the distribution is close to canonical (see below however), but also that 
its numerical value is very close to the one obtained from the model of Ref. 0- Indeed, the observed value of q in 
the present model was q = 0.1255 for TV = 8, which should be compared to q 0.1 obtained by BMWG for the same 
number of particles. 

In the following section we will be observing the temperature of the system via the average momentum square of 
the particle. Here I show that the mismatch between T$ and To to be found in that section can alternatively be 
attributed to the noncanonicity of f(p). Indeed, using the integral 

i/VI 1 ' 2of»+» r(j + 6+§) 



where T(u) is the gamma function and n — 0, 1,2, . . ., we obtain, after using the property T(u + 1) = uT(u) and 
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reading the values of a, b from Eq. I|34|) . 



/2E AT? 

</>=/ ^_/f{ P )d P =— (36) 



(compare with Eq. (|39|l below). Note that this result differs from the one obtained using the Boltzmann factor, 
namely (p 2 ) = 4E/3N (compare Eq. (|40[l ). This difference, though small for large N, is ultimately the motivation 
for BMWG to conclude that their corrected temperature is indeed an improvement over the "standard" one (Tq, 
obtained from Boltzmann's principle), which shows additionally that, contrary to the discussion that accompanies 
Fig. 6 of Ref. Q , the small departure from canonicity of the momentum distribution is a pertinent property of these 
small systems. 



B. Temperature 

Using the same weak coupling condition adopted in the last section, one can neglect the interaction Hamiltonian 
and compute both the phase space volume and area (this time for the whole system) through a procedure similar to 
the one leading to Eqs. (|32|l and (|33ll . One then finds, respectively, 

3(E) =5£( 3JV + 4 )/ 4 , (37) 
n(E) = dE 3N /\ (38) 

from which the temperatures T$ and To follow trivially, viz. 

* ~ V dE J ~ 3iV + 4 ~ ¥ 1 + 4/37V' (39) 

/.. (^m-) ; is if. (40) 



dE J 3JV 3 

where e = E/N is the total energy per degree of freedom and k was set to unity (this harmless simplification will be 
adopted throughout the remaining part of the paper). In what follows, the numerically observed temperature T b s is 
defined as the time-average v 2 , which should coincide with the ensemble average (v 2 ) [and hence with T$, cf. Eq. l|2"I)l] 
for ergodic systems. 

In Fig. [5]the observed temperatures from the numerical simulations with e = 1.0, g — 0.1 and various values of N 
are presented, as well as the analytic predictions Eqs. (|39|l and (|40|l . The agreement with T$ is evident, and we can 
see clearly from this figure that both temperatures should coincide in the thermodynamical limit 4] . In order to verify 
that these results can in fact explain the discrepancy found by BMWG between the "standard" form of Boltzmann's 
principle [which gives rise to Eq. (|19fl ] and their numerical simulations, the relative error £ — (T G ^ S — Tq^ /Tq was 
computed both for the data available here and in that reference (the temperatures from the latter are Tq w 14.95 
and T b s ~ 11.78), giving £a ~ —15.2% and £b ~ —21.2%, respectively. Note that not only the magnitude of these 
relative deviations are quite close, but also their sign is the same, i.e. the general trend of Tq is to overestimate 
the final equilibrium temperature, just as BMWG found (see also [l5|V The additional error £a — £b ~ 6% can 
be further explained from the fact that the phase space area in the approximate temperature expression of BMWG 
involves the degrees of freedom of the booster only, neglecting the contribution of x and p which tends to increase Tq 
and thus to decrease £ . This is the source of the "10%" deviation from the observed temperature mentioned before 
the conclusions of Ref. |5j ■ 



V. CONCLUSIONS 



It is important to emphasize that the purpose of the present contribution was not to determine whether the volume 
or the area entropy is the correct starting point of thermostatistics for small systems. Although the former bears an 
important thermodynamic property that is independent of the system size (known as "adiabatic invariance," cf. [l6| 
and references therein), these very small and isolated systems such as the ones studied by BMWG are usually far 
removed from experience, and we do not know whether our empirical and macroscopic thermodynamics acquires a 
different form in this case. Experiments with small clusters of atoms is a field that can potentially settle this question 
|17| , and a theoretical comparison between the two entropies above has already been carried out in this context |18| . 
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FIG. 2: Comparison between the observed temperatures of the particle T b 3 from the numerical simulations (full circles) 
and from the analytic predictions, Eqs. 1391 and 1401 (continuous and dashed lines, respectively). The observed equilibrium 
temperature T b 3 was obtained from the asymptotic value of the time-average v 2 after allowing a suitable relaxation and 
observation time (see text). 



Quite apart from elucidating the ad hoc definition of the BMWG temperature [cf. discussion following Eq. (|21|) ]. the 
present contribution has established that one does not have to resort to dynamical quantities such as time correlation 
functions in order to predict the equilibrium temperature of a small system. Such idea was put forward in the work 
of Ref. 5] in the form of a correction to the temperature that follows from Eq. for the particular case of a system 
composed of a particle coupled to a set of "irrelevant" degrees of freedom, leading its authors to conclude that a 
correction for the Boltzmann principle is necessary when one is dealing with small systems. The key identity that 
escaped the attention of BMWG is Eq. (|24[) . which establishes an exact relation between the equipartition temperature 
and the phase space volume, a quantity that can be computed much in the same way as the usual phase space area. 
This shows that, as long as we are working with the assumptions of Ref. the particular form of the Boltzmann 
principle in Eq. J2J) is sufficient for predicting the equipartition temperature, without any dynamical correction. 
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